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Resumen 


La temperatura es una variable trascendental en el cálculo de la 
evapotranspiración, el crecimiento, el desarrollo y el rendimiento de las 
plantas; en el estudio de la transmisión de plagas y enfermedades; en el 
pronóstico del clima; en la determinación del flujo de calor; en el cálculo 
de la presión real de vapor. Todos estos procesos son afectados por el 
calentamiento global. El objetivo de este trabajo es comparar los 
mejores resultados de dos modelos: uno de red neuronal artificial (RNA) 
backpropagation, y otro de series de Fourier. Se utilizaron datos diarios 
de temperaturas máximas (Tmax) y mínimas (Tmin) de las estaciones 
Santa Rosa 1, Ruiz Cortines, Batequis y Santa Rosa 2, del Distrito de 
Riego 075 Valle del Fuerte, Los Mochis, Sinaloa, México. En la RNA, 
1 484 vectores de datos fueron utilizados para entrenamiento, 
validación y prueba y 229 para pronóstico. Para el entrenamiento las 
variables de entrada de la RNA fueron: día juliano, longitud, latitud y 
altitud. Se obtuvieron 96 escenarios con una, dos y tres capas ocultas, 
con diversos números de neuronas en cada capa oculta. Con los 1 484 
datos, se obtuvieron los mejores ajustes para los modelos de series de 
Fourier para temperaturas máximas y mínimas, y se pronosticaron 229 
datos para las cuatro estaciones. Los mejores modelos de RNA 
backpropagation para el pronóstico de temperaturas máximas y 
mínimas diarias obtuvieron desempeños similares en comparación con 
los realizados por los mejores modelos de series de Fourier, para las 
estaciones de estudio. 


Palabras clave: Pronóstico, redes neuronales artificiales, temperatura 
máxima, temperatura mínima 


Abstract 
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Temperature ¡is a transcendental variable in aspects such as 
evapotranspiration calculation, growth, development and yield of plants, 
in the study of the transmission pests and diseases, in the weather 
forecast, in determination of heat fluxes, in the calculation of the real 
vapor pressure, all these processes affected by global warming. The 
objective of this work was to compare the best results of two models: 
one of artificial neural network (RNA) backpropagation, and another of 
Fourier series. Daily data of maximum temperatures (Tmax) and 
minimum (Tmin) of the Santa Rosa 1, Ruiz Cortines, Batequis and Santa 
Rosa 2 stations, of the Irrigation District 075 Valle del Fuerte, Los 
Mochis, Sinaloa, Mexico were used. In RNA, 1484 data vectors were 
used for training, validation and testing and 229 to forecasting. For 
training, the input variables of the RNA were: Julian day, longitude, 
latitude and altitude. Were obtained 96 scenarios with one, two and 
three hidden layers, with different numbers of neurons in each hidden 
layer. With the 1484 data, the best adjustments were obtained for the 
Fourier series models for maximum and minimum temperatures, and 
229 data were predicted for the four stations. The best RNA 
backpropagation models for the prediction of maximum and minimum 
daily temperatures obtained similar performances in comparison with 
those made by the best models of Fourier series, for the study stations. 


Keywords: Forecasting, artificial neural network, maximum 
temperature, minimum temperature 


Recibido: 11/08/2017 
Aceptado: 02/07/2018 


Introducción 


Las temperaturas máximas y mínimas son importantes en la estimación 
de la evapotranspiración de referencia (Hargreaves € Samani, 1985; 
Segura-Castruita 8  Ortiz-Solorio, 2017), para el cálculo del 
requerimiento de riego en los distritos de riego. También se emplean 
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para la detección local del cambio de la temperatura (Lobato-Sánchez 8 
Altamirano-Del-Carmen, 2017) que afecta la evapotranspiración de los 
cultivos. De igual manera presentan una marcada influencia en todos los 
procesos fisiológicos del crecimiento y desarrollo de las plantas 
(fotosíntesis o respiración), y son ¡importantes por los efectos 
perjudiciales en los cultivos (Elías € Castellví, 2001). 


Las redes neuronales artificiales son un modelo de pronóstico útil, dado 
la alta no linealidad que manejan. Cobaner, Citakoglu, Kisi 8 Haktanir 
(2014) para estimar las temperaturas mensuales mínimas, máximas y 
medias del aire en Turquía, se utilizaron una red neuronal artificial, un 
sistema de inferencia neuro-fuzzy y un modelo de regresión lineal 
múltiple o el planteamiento de Alexandridis 8 Zapranis (2013) de una 
red neuronal Wavelet para modelar la temperatura media diaria. 


Otros modelos también se han utilizado para estimar la temperatura, de 
forma horaria, diaria o mensual, por medio de regresiones lineales 
múltiples con variables predictivas como: latitud, longitud, altitud y 
continentalidad. Estos son planteamientos clásicos para la modelación 
de la temperatura (Monestiez, Courault, Allard, 8 Ruget, 2001 y 
Cristóbal, Ninyerola, 8 Pons, 2008). 


Por otra parte Thornley 8 France (2007) indican que la temperatura del 
aire tiene un comportamiento sinusoidal, por lo tanto es factible utilizar 
para su estimación un modelo de series de Fourier. El objetivo fue 
pronosticar las temperaturas máxima y mínima diaria del aire con redes 
neuronales artificiales para diferentes escenarios y funciones de 
transferencia (purelin y tansig) y comparar los mejores resultados con 
los obtenidos por el modelo de serie de Fourier para su uso en 
diferentes aplicaciones. 


Materiales y métodos 


La información diaria de temperaturas máximas y mínimas (abril 1997 a 
diciembre de 2001) fue obtenida de cuatro estaciones meteorológicas 
automáticas (EMAS), del Distrito de Riego 075 Valle del Fuerte, en Los 
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Mochis, Sinaloa, cuyos nombres (clave), latitudes, longitudes y altitudes 
son: Ruiz Cortínez (3843 II-2), 25% 39' 15”, 108% 45' 20”, 31 msnm; 
Batequis (3546 II-3), 25% 45" 49”, 108% 48' 41”, 32 msnm; AC Santa 
Rosa 1 (3765 III-1), 25% 45' 03”, 108% 57” 21”, 40 msnm; AC Santa 
Rosa 2 (9610 III-1) 25% 51' 16”, 108% 52' 03”, 61 msnm. 


Para obtener el modelo de serie de Fourier se utilizó el Curve Fitting Tool 
de Matlab (MathWorks Inc., 2001). El software para entrenamiento, 
validación, prueba y pronóstico de la red neuronal artificial fue el 
Toolbox de Matlab, nnet (Demuth € Beale, 2001). 


Modelo de serie de Fourier 


Con los datos de abril de 1997 a mayo de 2001 se generaron los 
modelos de serie de Fourier y de junio a diciembre de 2001 se utilizaron 
para la validación de estos. La variable de entrada fue día juliano, y la 
salida, temperatura máxima o mínima. 


La ecuación para la serie de Fourier presentada por Dubbel (1977) es: 
FO) = 40 +)Y;-, a¡cos(nwx) + b¡sen(nwx) (1) 


donde: aj, a; y b; son las constantes de la serie de Fourier, para / = 1 
oscilación u onda fundamental, y para ¡ = 2, 3,.....n, ondas superiores o 
armónicas superiores, w es la frecuencia fundamental de la señal, n es 
el número de términos (armónicas) en la serie, y en este caso: 1<n< 
8; x, es la variable dependiente. 


Redes neuronales artificiales 
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En la red neuronal artificial se aplicó el algoritmo backpropagation 
(Haykin, 1994). 


Para el entrenamiento, evaluación y validación se utilizaron los datos 
diarios de temperaturas máximas o mínimas de abril de 1997 a mayo de 
2001 y para el pronóstico de junio a diciembre del 2001. Las variables 
de entrada son: día juliano, longitud, latitud y altitud, y de salida: 
temperaturas máximas o mínimas. 


Las estructuras de los 96 diferentes escenarios para los entrenamientos 
de la red neuronal backpropagation se muestran en la Tabla 1, se 
combinó número de capas y número de neuronas en éstas, se utilizaron 
dos diferentes algoritmos de entrenamiento. 


Evaluación estadística de los resultados de los modelos 


Se utilizó el error estándar promedio (RMSE, por sus siglas en inglés), 


R a 
su ecuación es: RMSE = (Berti) , donde a, es el dato estimado por 


el modelo, t; es el dato observado y N es el número total de 
observaciones, (Cai, Liu, Lei, 8 Pereira, 2007). 


Tabla 1. Escenarios 1, 2, 3 y 4, estructura de construcción para éstos, 
con cambios en: función de transferencia (tansig, purelin), algoritmos 
de entrenamiento (trainbr, trainlm), número de capas ocultas, y número 
de neuronas en éstas 


Cuatro entradas para la RNA (día juliano, latitud, longitud y altitud) 


Escenarios 1 (el) Escenario 2 X(e2)> Escenario 3 (e3) Escenario 4 X(e4> 
Salida Tmax (función): Salida Tmax Salida Tmax Salida Tmax (función): 
f(tansig) para (el) (función): (función): f(purelin) f(purelin) 
f(tansig) para (e3) 
Salida Tmin (función): Salida Tmin Salida Tmin Salida Tmin (función): 
f(tansig) para (el) (función): (función): f(purelin) f(purelin) 
f(tansig) para (e3) 
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Algoritmo de Algoritmo de Algoritmo de Algoritmo de 
entrenamiento entrenamiento entrenamiento entrenamiento trainbr 
trainlm trainlm trainlm 
£neuronas» £neuronas» £neuronas)> £neuronas)> 
1 capa oculta 1 capa oculta 1 capa oculta 
fítansig) j fítansig) fítansig) 
LA Si (4 4 
(85 e (85 (85 
127 - (12) (12) 
(16) Ñ (16) (16) 
(20) 5 (20) (20) 
1247 E (24) (24) 
(28) - (28) (28) 
(32) - (32) (32) 
(36) = (36) (36) 
(407 - (40) (40) 

2 capas ocultas, 2 capas ocultas 2 capas ocultas, 2 capas ocultas 
fítansig)xKtansig)y Ktansigyxftansig) fítansig)>xfftansig)y fítansig)xfítansig) 
(4x4) (8x4) (4x4) (8x4) 
(8x8) (12x4)> (8x8) £12x4) 
(12x12) (16x8)> (£12x12) (£16x8) 
(16x16) £20x8)> (£16x16) £20x8) 
£20x20y £24x16) £20x20) (£24x16) 
(24x24) (28x16) £24x24) (£28x16) 
£28x28y (£32x16) (£28x28) (£32x16) 
£32x32y £36x20)> (£32x32) £36x20) 
(36x36)$ £40x20)> (£36x36) £40x20) 

£40x40)> - £40x403 - 

3 ocultas 2 capas ocultas, 3 capas ocultas, 
fítansig)xfítansigyxf 7 Ktansig xK tansig y fítansig»xKtansigyxKX 
£tansig) tansig)> 
(4x4x4> - 8x4 (8x8x4) 
(8x8x8> - 12x4 (£8x12x4) 

£12x12x12) - 16x8 (£8x16x4) 
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(£16x16x16) - 20x8 £8x20x4) 
£20x20x20) S 24x16 £8x24x8) 
£24x24x24) - 28x16 (£8x28x8) 
(£28x28x28) - 32x16 (£16x32x16) 
(£32x32x32) - 36x20 (16x36x16) 
(36x36x36) a 40x20 (£16x40x16) 
£40x40x40> - - - 


Tmax: temperaturas máximas; Tinin: temperaturas mínimas; trainlm: Algoritmo de entrenamiento Levenberg- 
Marquardt; trainbr: Algoritmo de entrenamiento de regularización Bayesiana; tansig: Función de 
transferencia tangente-sigmoidea; purelin: Función de transferencia de limite rígido. 


Resultados y discusión 


Entrenamiento, evaluación y validación (1 484 datos), 
para los 96 escenarios de la RNA 


Se determinaron todos los escenarios (Tabla 1), los rangos de variación 
de los RMSE de estos, para la variable temperatura máxima fueron de: 
1.2540 (e4) (20) a 2.4858 (e4) (40x20); 2.1790 (e3) (40x40) a 
2.2765 (e4) (8x12x4); 2.3280 (el) (40x40x40) a 2.4277 (e4) 
£16x40x16) y 2.5093 (e3) (40) a 2.6276 (e4) (16x32x16)?, para las 
estaciones Santa Rosa I, Ruiz Cortínez, Batequis y Santa Rosa 2, 
respectivamente. 


Para la temperatura mínima fueron de: 2.2372 (el) (£36x36) a 2.2554 
(e4) (32x16); 2.0152 (e3) (40x40) a 2.17 (e4) (16x40x16); 2.1523 
(e3) (36x36) a 2.2934 (e4) (8x8x4) y 2.01 (e3) (40x20) a 2.0809 
Lfe4) £32x16), para las estaciones Santa Rosa 1, Ruíz Cortínez, 
Batequis y Santa Rosa 2 en esa disposición. 
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Pronóstico (229 datos) con la RNA 


Los rangos de los RMSE de todos los escenarios (Tabla 1), primero para 
la temperatura máxima, luego para la temperatura mínima, y para las 
estaciones: Santa Rosa 1, Ruiz Cortines, Batequis y Santa Rosa 2, en 
ese orden fueron de: 2.4448 (e4) (12x4) a 2.5522 (el) (4x4x4); 
2.4697 (e3) (4x4) a 2.5420 (e2) (8x4); 2.4915 (e1) (4) a 2.5246 
[el) (4x4x4); 2.8649 (e1) (12) a 2.9574 Xe1) (4x4), y 2.0223 (e4) 
£16x32x16) a 2.1713 (e2) (36x20); 1.9922 (e4) (16x32x16) a 2.09 
Le4) (24x16); 2.0346 (e4) (8) a 2.1776 (e2) (12x4); 2.3035 (el)> 
(8) a 2.3865 (e3) (32x16). 


De todos los escenarios presentados para el entrenamiento y el 
pronóstico de las temperaturas máximas y mínimas diarias son mejores 
estimadores los escenarios (e3) para el primero, (e4)3 y Xel) para el 
segundo, para estos se aplica en la capa de salida la función de límite 
rígido (purelin), a excepción de (ely, y en las capas intermedias la 
función tansig. Sahin (2012) realizó estimaciones de temperaturas 
medias mensuales del aire con sensores remotos y redes neuronales 
artificiales y obtuvo RMSE de 1.254 y 1.263 K en sus mejores ajustes y 
utilizó las mismas funciones de transferencia, en la salida y en la capa 
oculta. La mayoría de los RMSE en el presente estudio no son menores 
de 2*C, ya que se utilizaron valores de temperatura diaria. 


Las diferencias de los RMSE de todos los rangos de los escenarios 
presentados tanto para el entrenamiento como el pronóstico no son 
mayores a 0.16%C para ambas temperaturas (máxima y mínima) y en 
las cuatro estaciones (Tabla 2), por lo que se puede utilizar cualquier 
escenario. Es la excepción la temperatura máxima de la estación Santa 
Rosa 1 (1.2%C). 


Tabla 2. Diferencias de los rangos de variación de los RMSE obtenidos 
para las variables Tmax Y Tmin- 


Entrenamiento (1484 
Estación datos) 


Tmax Tmin T max T min 


Pronóstico (229 datos) 
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Santa Rosa 1 1.2318 0.0582 0.1074 0.149 

Ruiz Cortines 0.0975 0.1548 0.0723 0.0978 
Batequis 0.0997 0.1411 0.0331 0.143 

Santa Rosa 2 0.1183 0.0709 0.0925 0.083 


Tmax: temperaturas máximas; Tinin: temperaturas mínimas. 


Los mejores escenarios de los entrenamientos no coinciden con los 
mejores escenarios de los pronósticos (Tabla 3). Además en la mayoría 
de los mejores entrenamientos, para ambas temperaturas, el mayor 
número de neuronas en las capas ocultas tienen un RMSE cercano a 
cero, pero no necesariamente los entrenamientos con el mayor número 
de capas ocultas serán siempre los que obtengan los mejores ajustes, 
solo uno de los escenarios con tres capas ocultas presentó uno de los 
mejores ajustes (estación Batequis). También se observa que a 
diferencia de los entrenamientos, en los pronósticos la mayoría de los 
mejores ajustes se presentan cuando el número de neuronas en las 
capas ocultas son las menores, esto corrobora lo que comentan Demuth 
8: Beale (2001) y Tymvios, Michaelides 8 Skouteli. (2008), que uno de 
los problemas que ocurre cuando se está entrenando una RNA es que se 
sobreajusta el conjunto de entrenamiento y no generaliza bien para un 
nuevo conjunto de datos (pronóstico), como sucede en el presente 
trabajo. 


Además (Tabla 3) no hay un escenario determinado que haya sido el 
mejor, tanto en el entrenamiento como en el pronóstico. 


Lo que interesa son los mejores escenarios del pronóstico, ya que estos 
garantizan que los datos estimados estén más cercanos a los 
observados (RMSE cercanos a cero). Como los mejores escenarios del 
pronóstico no se corresponden con los mejores del entrenamiento, se 
tomó la decisión de trabajar con los del pronóstico que se presentan en 
las dos últimas columnas (Tabla 3). 


Tabla 3. Mejores escenarios con sus respectivas neuronas en la(s) 
capa(s) oculta(s) y RMSE, para el entrenamiento y pronóstico. 


a Mejores escenarios con su error 
Estación 


Entrenamiento Pronóstico 
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T max T min T max Tmin 
(e4) (e4) 
El die da [e1) £36x36), 2.24 (12x4), (16x32x16), 
: 2.44 2.02 
Rúl (e3) [e3) (40x40>, 2.02 (e3) (e4> 
Cortinés (4x40)>, (4x4),  (16x32x16), 
2.18 2.47 1.99 
(el) (el) 
Batequis (40x40x40), (e3) (36x36), 2.15 44), O 
2.33 2.49 : 


(eL, 
Santa Le3) (40), (el) (8), 
o e [e3) 40x20), 2.01 az), San 


Tmax: temperaturas máximas; Tmin: temperaturas mínimas 


Comparación de los mejores escenarios del pronóstico 
de RNA con los modelos de series de Fourier de mejor 
ajuste 


Los mejores escenarios del pronóstico de las RNA se compararon con los 
mejores modelos de series de Fourier obtenidos para cada una de las 
estaciones, en la Figura 1 (Santa Rosa 1) se muestra el comportamiento 
de las diferencias de los valores observados de temperatura máxima con 
los valores pronosticados con cada uno de estos dos modelos. Estas 
diferencias presentan las mismas tendencias y están alejadas de la línea 
del valor cero en la misma proporción. Además se observa que el rango 
de las diferencias está entre 10% C y -6% C. Esto que se comentó se 
presenta también en las otras tres estaciones y para la temperatura 
mínima (Figura 2). 
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101 E 

al | 


-RNA[e4X12x4)-datos observados, RMSE = 2.4448 | 
= =f(x)-datos observados, RMSE = 2.6213 ij y 


Diferencias temperatura máxima 


En 00 a - E ca 
atos de pronóstico 
Figura 1. Diferencias de los modelos de RNA y serie de Fourier con los 
datos observados de temperaturas máximas para Santa Rosa 1 AC; 
modelo de Fourier con n = 5: f(x) = 31.86 -— 0.6201 cos(0.01727x) + 
5.871 sen (0.01727x) + 0.8935 cos(2*0.01727x) - 0.0709 
sen(2*0.01727x) - 0.745 cos(3*0.01727x) + 0.6926 sen(3*0.01727x) 
+ 0.3822 cos(4*0.01727x) + 0.1923 sen(4*0.01727x) - 0.227 
cos(5*0.01727x) + 0.1547 sen (5*0.01727x). 
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Figura 2. Diferencias de los datos observados de temperaturas mínimas 
con los modelos de RNA y serie de Fourier, para Ruiz Cortínez; modelo 
de Fourier con n = 8: f(x) = 15.98 - 2.83 cos(0.0172x) + 8.476 sen 
(0.0172x) - 0.5298 cos(2*0.0172x) - 1.121 sen(2*0.0172x) - 1.142 
cos(3*0.0172x) + 0.3966 sen(3*0.0172x) + 0.1746 cos(4*0.0172x) - 

0.6291 sen(4*0.0172x) + 0.1938 cos(5*0.0172x) - 0.1892 sen 
(5*0.0172x) - 0.04237 cos(6*0.0172x) + 0.09564 sen(6*0.0172x) + 
0.2258 cos(7*0.0172x) + 0.05545 sen(7%*0.0172x) - 0.4136 
cos(8*0.0172x) - 0.02458 sen(8*0.0172x). 


En la estimación de las variables temperatura máxima y mínima, los 
mejores modelos de series de Fourier solo para dos estaciones, están 
descritos debajo de las Figuras 1 y 2. Los RMSE obtenidos, para el 
ajuste (1484) y la validación (pronóstico) del modelo (229), para la 
temperatura máxima fueron: Santa Rosa 1, con n=5, 2.553, 2.6213; 
Ruiz Cortines, con n =7, 2.266, 2.5321; Batequis, con n=7, 2.417, 
2.6847 y Santa Rosa 2, con n =7, 2.555, 3.0855. Y para la temperatura 
mínima: Santa Rosa 1, n = 7, 2.379, 2.1415; Ruiz Cortines, n = 8, 
2.163, 2.1759; Batequis, n = 8, 2.289, 2.1822 y Santa Rosa 2, n=7, 
2.135, 2.5005. 


Al comparar los RMSE (Tabla 3) de las RNA con los del modelo de series 
de Fourier para temperatura máxima en las cuatro estaciones para el 
pronóstico, los primeros resultan mejores, por centésimas, lo mismo 
sucede al comparar los RMSE para la temperatura mínima. Al tener en 
cuenta que las temperaturas se registran hasta decimas de grado, las 
RNA presentan un mejor desempeño (menores RMSE) en máximo tres 
décimas, que el modelo de Fourier para todas las estaciones. No 
obstante que para entrenar este tipo de red neuronal (backpropagation) 
se necesitan conocer cuatro variables de entrada (día juliano, latitud, 
longitud y altitud) y para el modelo de series de Fourier únicamente la 
serie de tiempo en días. 


Lo obtenido en este trabajo corrobora lo encontrado por Sahin (2012), 
que propone a la red neuronal artificial como un método alternativo para 
estimar la temperatura del aire con suficiente exactitud. La importancia 
de los modelos de series de Fourier obtenidos es como lo comentan 
Coronas €: Baldasano (1984), que las ecuaciones son representativas de 
forma anual para donde se generaron y su ventaja, es que representan 
una gran cantidad de datos en una sola ecuación y se incorporan nuevos 
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datos fácilmente. Finalmente, Almonacid, Pérez-Higueras, Rodrigo € 
Hontoria (2013) también realizaron una comparación de tres métodos 
para estimar temperaturas horarias para localidades españolas: el de 
redes neuronales artificiales proporcionó los mejores resultados, la 
función de transferencia que utilizaron fue purelin en la capa de salida y 
tangente sigmoidea en la capa oculta, obteniendo RMSE en un rango de 
0.53 - 1.20% C, estas funciones fueron las que mejor desempeño 
obtuvieron en el presente trabajo tanto en los ajustes, como en los 
pronósticos. 


Conclusiones 


En el pronóstico de las temperaturas máximas o mínimas en el distrito 
de riego 075 para cualquiera de sus aplicaciones, como el cálculo de la 
evapotranspiración de referencia, se pueden utilizar cualquiera de los 
escenarios de las RNA. Las funciones de transferencia que ofrecen 
mejores desempeños en el ajuste para el entrenamiento y/o pronóstico 
de la temperatura máxima o mínima diaria con la RNA backpropagation 
fueron: tangente sigmoidea en la capa oculta y purelin en la capa de 
salida. Los mejores modelos de RNA para el pronóstico de temperaturas 
máximas y mínimas diarias obtuvieron desempeños similares en 
comparación con los realizados por los mejores modelos de series de 
Fourier, para las estaciones de estudio. 
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